Monte Carlo simulation based on dynamic disorder model in organic semiconductors: 

Prom bandlike to hopping transport 
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The dynamic disorder model for charge carrier transport in organic semiconductors has been ex- 
tensively studied in recent years. Although it is successful on determining the value of bandlike 
mobility in the organic crystalline materials, the incoherent hopping, the typical transport charac- 
teristic in organic semiconductors, cannot be described. In this work, the decoherence process is 
taken into account via a phenomenological parameter, say decoherence time, and the projective and 
Monte Carlo method is applied for this model to determine the waiting time and thus the diffusion 
coefficient. We find the type of transport changes from bandlike to incoherent hopping with a sufB- 
ciently short decoherence time, which indicates the essential role of decoherence time in determining 
the type of transport in organics. We have also discussed the spatial extent of carriers for difi'erent 
decoherence time, and the transition from delocalization (carrier resides in about 10 molecules) to 
localization is observed. Based on the experimental results of spatial extent, we estimate the deco- 
herence time in pentacene has the order of Ips. Furthermore, the dependence of diff'usion coefficient 
on decoherence time is also investigated, and corresponding experiments are discussed. 

PACS numbers: 72.80.Le, 72.20.Ee, 03.65.Yz, 63.50.-x 
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I. INTRODUCTION 

It has been long-termly known that, the two differ- 
ent types of charge carrier transport, namely, incoher- 
ent hopping and bandlike tunneling, coexist in organic 
semiconductors, and both of them play an essential role 
simultaneouslyfii^ At the very beginning, when the con- 
ducting conjugated polymer was firstly found in 1970s, 
the explanation for the underlying mechanism of the elec- 
tric conductivity was based on the picture of solitons and 
polaronsj^ It is the strong electron-phonon (e-p) coupling 
in organic materials that makes the electrons or holes self- 
trap in a lattice distortion. This charged polaron is the 
main carrier, which is initially realized to be completely 
localized in an individual molecule. In this context, it 
was smoothly concluded that, the character of charge 
transport between organic molecules should be incoher- 
ent hopping of polarons. However, in the last decade, 
plenty of experiments reported that there is a region that 
bandlike tunneling works. ^'t^— For example, there exists 
an abnormal region around 100-200K that the mobility 
decreases as temperature increases (i.e., negative temper- 
ature coefficient of mobility). This indicates the type of 
carrier transport in this region is bandlike and coherent,- 
as the charge is found to be delocalized in a region of 
several molecules (about 10 molecules from the experi- 
ment of electron spin resonance^) rather than a single 
molecule. This result is somehow contradiction to the 
traditional understanding of localized polarons. In ad- 
dition, it was also stated by a recent experiment on the 
ultrafast initial carrier dynamics that, the intramolecular 
microscopic dynamics should be significant in the charge 
transport/^ Up to now, to describe the dualistic coherent 
and incoherent transport is still an open subject. 

To understand the new experimental findings, a num- 
ber of theoretical works based on the microscopic mod- 



els, such as dynamic disorder model,™"— Su-Schrieffer- 
Heeger (SSH) model^^ and Anderson model?^ were 
addressed to account for the coherent and incoherent 
transport)^""— and the coherent length and diffusion 
constan1i2£ have also been studied. Based on Holstein- 
Peierls model, Troisi and his coworkers highly com- 
mended the mechanism of dynamic disorder. They pro- 
posed that, due to the scattering with phonons, the 
carrier will be localized into a small region, and thus 
the mobility decreases. Their theory may be applied 
to the temperature dependence of mobility in organic 
crystal but has not comprehensively described that in 
organic semiconductors Meanwhile, under sufficiently 
long time evolution, the model seems failing to capture 
the basic bandlike characteristic.^^ More importantly, the 
Ehrenfest method used by Troisi has some fundamental 
drawbacks, and the typical one is that the superposi- 
tion principle is deviated.-- This is because the deco- 
herence process has not been appropriately taken into 
consideration^ while for the simulation of incoherent 
hopping transport in organic semiconductors, the deco- 
herence process should be of actual importance. 

In the common sense, the role of phonons is twofold: 
One is to provide the energy to help the carrier hop- 
ping (hopping transport) and the other is to scatter with 
carriers to obstruct them (bandlike transport). These 
two mechanisms might coexist in organic semiconduc- 
tors, which makes the discussion on this issue very com- 
plicated. In this work, we apply the decoherence process 
into Troisi's dynamic disorder model^^ and then embed it 
into the Monte Carlo simulation to quantitatively evalu- 
ate the diffusive coefficient within organic molecules. Our 
aim is to provide a generalized and efficient way to merge 
the above two roles of phonons into a single framework 
and find the transition between bandlike and incoherent 
hopping transport. The paper is organized as follows. 
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The next section will be the introduction of the present 
theory. The main results will be in the third section, and 
the final section is for the conclusions. 



picture of incoherent hopping should be addressed, which 
is one of the main subject in this work. 



II. INTERPLAY OF INCOHERENT HOPPING 
AND COHERENT EVOLUTION 

In this section, we will introduce the present theory. 
As the basis, incoherent hopping and thus the decoher- 
ence process during intermolecular hopping will be firstly 
described by introducing a phenomenological parameter, 
i.e., decoherence time. Then the details of the calculat- 
ing procedure of diffusion coefficient will be described. 
At last, the dynamic disorder model will be introduced 
into the theory. 



A. Incoherent hopping 

We first discuss the incoherent hopping, which is the 
motivation we consider the process of decoherence in 
this work. Originally, in order to study the incoherent 
hopping, Bassler applied the kinetic Monte Carlo sim- 
ulation with the so-called Miller-Abrahams (M-A) for- 
mula embeddedi^ Based on the Gaussian disorder model 
(GDM), he was able to evaluate the relationship be- 
tween the mobility of carriers and temperature, elec- 
tric field, and energy disorder. Especially, he found 
for the temperature dependence of mobility a scaling 
^ exp(— r^^), which is quite different from the tra- 
ditional understanding..23. Following him, there are lots 
of work investigating the infiuence of carrier density^^i 
electron-electron interactions and so oni^ Meanwhile, 
the other theories, such as the molecular dynamicSfiS, the 
percolation theoryj.2L22, and those considering the trap- 
ping mechanism ;^^!^" are applied extensively to comple- 
ment the discussions of this incoherent hopping. Espe- 
cially, the dynamic disorder model has already been put 
into the Monte Carlo simulation under the displaced har- 
monic oscillator approximation,^ which is somehow dif- 
ferent from the present treatment. 

Fundamentally, the M-A formula, coinciding with the 
detailed balance, is related to a single phonon process, 
such that the energy difference between initial and fi- 
nal site must be comparable with the highest energy of 
phonons that could effectively interact with carriers 
Hence, the hopping mechanism based on M-A formula 
could be described like this:-- initially a few of vibra- 
tional modes (phonons) are excited, such that the carri- 
ers could be heated up to make the energy levels of two 
neighboring sites coincide, and then the phonon transit 
its energy to the carrier to assist it hopping. This mech- 
anism correspondingly provides an intuitive picture to 
understand the application of M-A formula in the case of 
low carrier's density^^^ but in common cases, the vibra- 
tional modes are quite diverse, which means the single 
phonon process loses its efficiency. Hence, a converted 



B. Process of decoherence 

In organic materials, charge carriers move in a com- 
plicated environment, which is mainly composed of the 
disorders of molecular configurations and vibrations, and 
decoherence happens incessantly. Normally, the high fre- 
quency modes of intermolecular vibration (fast interac- 
tion) could be realized to be the source of decoherence, 
while the lower ones (slow interaction) is referred to a 
thermal reservoir as in dynamic disorder model discussed 
in belowJ^ The typical decoherence time is of the or- 
der Ips in organic materials with good conductivity;^ 
such as in C60, where the frequency of center-of-mass 
motion is observed to be 1.2THz,'^- but in more amor- 
phous materials, the value should be smaller. The com- 
mon method treating decoherence in molecular systems 
is to add some fiuctuation and dissipation to the model 
and evaluate the corresponding Green's function. ■^^ Ob- 
viously, this method is invalid for the dynamic disorder 
model, in which the motion of lattice sites must be cal- 
culated within real spacetime. So here we will describe 
a more efficient picture for the decoherence process^^ 
which is straightforward from a so-called coarse graining 
methodi^ 

To this end, we take two molecules indexed by 1 and 
2 as an example, as shown in Fig. [1] Initially, there is a 
carrier residing in molecule 1 as a localized wavepacket. 
It will spread following the equation of motion. 



(1) 



where we have used the tool of density matrix 
Pii Pv2 \ j-^^, convenience to understand deco- 
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herence. 7? in ([ij is the Hamiltonian considering the 
simplest hopping term of polarons, namely. 



77 = ri2(CjC2 + /i.c.), 



(2) 



where C\ {Ci) creates (annihilates) a polaron at i(= 1, 2)- 
th molecule, t\i the intermolecular overlap integral of 
electronic wavefunctions. 

When the frequency of vibration is high enough, the 
vibrational modes make the decoherence happen very 
quickly, •^^'•^^ such that we can define a decoherence time 
td, within which the carrier will lose its coherence and 
the off-diagonal terms pn and p2\ carrying the phase in- 
formation of the system vanishfSi Normally, td is of the 
order of 1/ti_2'^ After i^, the pure density matrix p be- 
comes a statistical mixture, and due to the vanishing of 
overlap between 1 and 2, the system could be realized to 
be a classical ensemble with pn and pii the probability 
that the carrier is residing in molecule 1 and 2, respec- 
tively. One can then produce a random number ?' e [0,1) 
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FIG. 1: Schematic for the process of incoherent hopping. Ini- 
tially, there is a polaron residing in molecule 1. After time td, 
its wavepacket will spread to molecule 2. Then the correlation 
between the two molecules is quenched due to the decoher- 
ence, and it will be decided whether one try of the incoherent 
hopping is successful. If yes, the polaron will hop to molecule 
2, otherwise the whole process restarts. 

as the procedure in Monte Carlo simulation. If r < p22 , 
the carrier will be chosen to reside in molecule 2 after td, 
and one step of the incoherent hopping finishes. Other- 
wise, the carrier remains at molecule 1, and the above 
process restarts. This provides the generalized picture of 
decoherence. 



C. Evaluation of waiting time 

In the common studies of charge transport in organ- 
ics, the coherent motion of carriers could be simulated 
within both adiabatic^^ and nonadiabati o^^' dynami- 
cal method. As well, when one wants to study the inco- 
herent hopping only, the Marcus theory based on lattice 
relaxation and level crossing is frequently used4£ Green's 
function was also evaluated to distinguish incoherent and 
coherent processJ^ In this subsection, we will develop 
such a method that merges both coherent and incoher- 
ent process on the basis of decoherence as we described 
above. The central task is that, when a carrier is initially 
residing in a molecule, we need to find a way to deter- 
mine the probability that a carrier hops out of the initial 
molecule. 

The procedure of our treatment is as follows. Initially, 
there is a carrier localized at molecule 1 (or the central 
site of the chain, as discussed below), that is, the initial 
state of the system is 

\I) = \l)^Cl\0). (3) 

The wavefunction of this carrier will then spread follow- 
ing Eq. ([1]). Since the carrier loses its coherence every 



td as discussed above, it will choose a molecule (e.g., 
molecule 2) to reside in based on the probability at this 
molecule at td- It means that, we calculate the probabil- 
ity of carrier at molecule 2, i.e., P2{td){= {td\C\C2\td)) 
with \td) the state of system at time td, and then make 
P2{td) the hopping probability that the carrier hops out 
of the initial molecule. This treatment is almost the same 
with the process of decoherence as we discussed above. 

Obviously, when td is small enough, it is extremely dif- 
ficult for a carrier to hop out of the initial molecule by 
one try. So the key step of the present theory is that, 
when the try fails, namely the carrier keeps at the initial 
molecule, the wavefunction of system will be projected 
onto the initial molecule, and continue evolving. That is, 
a new round of evolution starts with the renewed projec- 
tive initial state 

|/„e») = |l)(l|id). (4) 

Again we calculate P2{2td) after normalizing \Inew) by 
making the above steps recur. Of course, if necessary, td 
could also be randomly selected. Now, we have in hand a 
series of time-dependent P2{ktd), which is different from 
those time-independent hopping rates in GDM, so the re- 
maining question is to determine the waiting time t^ for 
a carrier residing in the initial molecule. Following the 
usual idea of Monte Carlo simulation, we sum up P2{ktd) 
and make the hopping happen when the summation re- 
sult exceeds a given random number, i.e., 

£P2(H,) >-lne, (5) 

k=l 

where ^ is a random number uniformly distributed be- 
tween and 1. In case for simplicity, the right hand side 
of ([5]) could also be replaced by 1. Then the waiting time 
is 

tw — kyjtd- (6) 

The above procedure could obviously be generaUzed for 
different systems beyond two molecules. In the following, 
we will apply it to the dynamic disorder model. 

Before moving on to the next subsection, we would like 
to say that, the projective method we describe here is al- 
most the same with the coarse graining method, with the 
only difference that the projection is acting within the 
real space other than energy space.-^^ So that, the basic 
drawback of Ehrcnfcst method we will use shortly, such 
as the deviation of superposition principle, is overcome in 
some sense. Meanwhile, the sudden switching method of 
electronic states in molecular dynamics is used for refer- 
ence of the Monte Carlo simulation. Equivalently, as we 
will show later, one can also use a rate equation with the 
off-diagonal decay term to calculate the diffusion process, 
which will give similar results with the present theory in 
some case. But for a multi-site system, since the veloc- 
ities of decay of the off-diagonal term are not exactly 
the same, the rate equation becomes quite complicated 
and inefficient in practise. This is why we consider the 
projective procedure here. 
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D. Dynamic disorder model 

The dynamic disorder model recommended by Troisi 
is a quite simple but efficient one in the study of charge 
transport in organicsiiSiii The basic idea is to use a one- 
dimensional SSH-like Hamiltonian to describe the cou- 
pling between charge carriers and intermolecular vibra- 
tional mode. Since the initial distribution of vibrational 
mode is dependent on temperature, the disorder intro- 
duced by this electron-phonon coupling is changeable 
with temperature, which implies the original meaning of 
the words " dynamics disorder" . The model Hamiltonian 
writes, 

H = Hcic + Hvib- (7) 

The electronic part is 

Hele = - Xlt^ ~ a{Uj + i - Uj)]{c]^^Cj + h.c), (8) 

j 

where c|(cj) creates (annihilates) a carrier on the j-th 
site. Uj represents the displacement of the j-th site, r 
is the transfer integral and a is the electron-lattice cou- 
pling constant. The vibrational part of Hamiltonian ([7]) 
is described as 

^^vib = y ~ )' + Y I] 

j j 

where K is the elastic constant between neighbor sites 
and M the mass of a site. The parameters are chosen to 
the similar ones as in [10]. That is, K = 14500amu ps~^, 
M = 250amu, a = 100 - 995cm- VA, t = 30-300cm-\ 
and the lattice constant is 4A. The number of site for 
the chain will be set to sufficiently large, such that the 
diffusive wavepacket of carrier can not reach the end of 
the chain. 

Based on this Hamiltonian, we could use the Ehrenfest 
method to calculate the diffusion of an initially localized 
wavepacket of carrier The initial state could have sev- 
eral forms and result in slightly different trajectories, but 
the basic behavior does not change4i In this work, the 
initial electronic state we choose is to locate a carrier at a 
single site, e.g. the central site of the chain. And the ini- 
tial condition of vibrations is under thermal equilibrium, 
i.e., {ui} and {iii} will be randomly chosen in a Gaus- 
sian distribution with variance ksT / K and ksT /M , re- 
spectively, with T the temperature. Then we apply the 
time-dependent Schrodinger equation for electronic part 
and Newtonian equation for vibrational part to compute 
the time evolution of the whole system. 

In the original treatment of Troisi, the evolution time 
should be long enough to ensure the saturation of diffu- 
sion coefficient. But in the present work, we calculate the 
evolution for each round to the time td and then make 
a decision whether the carrier could hop out of the ini- 
tial site. The whole calculating procedure for electronic 



part is almost the same as described in above subsec- 
tions. However, as we are now facing a multi-site system, 
we treat the initial site as molecule 1 and the others as 
molecule 2, and after the hopping is determined, we need 
to produce another random number to determine which 
site the electron should go. On the other hand, the evo- 
lution of vibrational part will be successive during the 
decision of the final site of electron. The quantities we 
compute are the diffusion coefhcient defined as 

D = (r2)/2t^, (10) 

with r the distance between the initial site and the final 
site that the carrier hops to, whose mean value is the 
diffusion length as we will show in below. Of course, 
if necessary, the mobility could be calculated by Ein- 
stein relatioufiS but Einstein relation is violated under 
nonequilibrium condition;^ such that the validity to use 
it here should be doubted. So we did not plan to plot it 
in the present work. 

III. RESULTS AND DISCUSSIONS 
A. Comparison with no decoherence case 

In Fig. [21 we show the comparison between the present 
theory and others. Firstly, under the parameters that 
Troisi has used,— we change the value of td and cal- 
culate the same quantity with Troisi, i.e., the averaged 
squared displacement. The red lines are from our projec- 
tive method, while the blue ones are from the calculation 
based on the master equation. For the latter case, we 
add an exponential decay term with variance td for all 
the off-diagonal elements of the density matrix and cal- 
culate the evolution and the physical quantities we want. 
This treatment is an intuitive way to treat the decoher- 
ence at least in two-site systems. 

As shown in left panels of Fig. [21 following the td de- 
creasing, the diffusion becomes much slower than Troisi's 
results, say no decoherence case. Meanwhile, when the t 
and a reduce by one order, the reduction of diffusion co- 
efficient becomes quite large. These results clearly show 
that, when the decoherence is presence, the ability of 
charge transport becomes poorer, and as we will show 
shortly, the type transits from bandlike to hopping trans- 
port. On the other hand, the results from the present 
projective method are very close to those from the mas- 
ter equation, but the computational efficiency improves 
a lot. 



B. Transition from bandlike to hopping transport 

The most important characteristic difference between 
bandlike and hopping transport is that, for the former 
one the diffusion coefficient decreases with increasing 
temperature due to the thermal induced dynamic dis- 
order (say, a power law scaling ~ T"™), while for the 
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FIG. 2: Averaged squared displacement versus time for differ- 
ent parameters with (red and blue) and without decoherence 
(black). The red lines are from our projective method, while 
the blue ones are from the calculation based on the master 
equation. 
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latter case the diffusion coefficient increases witfi increas- 
ing temperature due to the thermal assisted mechanism 
cxp(— T^™)). In Troisi's original studies, only band- 
like transport was found for all the range of parameters. 
Here in Fig. [3l we show for different parameters the diffu- 
sion coefficient. Our results clearly show that, when r is 
small (r — 30cm~^, which is the typical value in organic 
small molecules), there is only positive temperature de- 
pendence of diffusion coefhcient. While when t is large 
(r = 300cm~^, which is the typical value in pentacene), 
there is only a negative and power law (m ~ 0.5) temper- 
ature dependence. In between (r = lOOcm"-'^), both pos- 
itive and negative temperature dependence emerge. On 
the other hand, as shown in Fig. El^b), when we change td, 
we find a critical region of the transition between bandlike 
(low temperature) to hopping (high temperature) trans- 
port, and the transition point is around lOOK. This is the 
most important finding in this work. The explanation of 
these results is that, for large r, the case is the same 
with that discussed in dynamic disorder modeljiS, while 
for small r, the carrier has little time to diffuse to other 
sites, but some phonon induced potential well could assist 
carrier to move. So in the latter case, temperature plays 
a positive role, and the type of transport becomes inco- 
herent hopping. More interpretation will be addressed 
shortly when we study the diffusion length. In addition, 
to get the sense of the influence of electron-phonon cou- 
pling, we also show the results with a = lOOcm^^/A in 
Fig. ^c). Based upon our calculation, the transition re- 
gion from bandlike to incoherent hopping is insensitive 
to a except some quantitative change, which means the 
dynamic disorder itself does not influence the transition 
of the two type of transport. 

As we have mentioned, the temperature dependence 
is mostly negative in crystalline organic semiconductors, 
whose mechanism should be bandlike transport. Here, 
applying the physics of decoherence, we obtain both pos- 
itive and negative temperature dependence for different r 
and temperature. For the typical value of r in pentacene, 
namely r — 300cm~^, it shows negative temperature de- 
pendence, which is the same with the result in dynamic 
disorder modelJ^ It implies that, in this material, band- 
like mechanism dominates the charge transport, and the 
role of decoherence is just to decrease the diffusion coef- 
ficient as we will discuss later on. However, the value of 
r in other organic materials, such as Alq3, is about two 
order smaller than that in pentacene, i.e., < BOcm^^, so 
the temperature dependence becomes positive. Hence, 
in common organic semiconductors under room temper- 
ature, the characteristic of hopping transport should be 
more frequently observed. 



C. Spatial extent of carriers 



FIG. 3: Diffusion coefficient versus temperature for different 
r, td, and a. The values of a is in the unit of cm^^/A. 



The spatial extent of carriers is a controversial but cru- 
cial issue in this field. Understanding of it is essentially 
helpful to clarify the interplay between bandlike and hop- 
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FIG. 4: Diffusion length (in the unit of the number of site) 
versus temperature for different td- t = lOOcm"^ and a = 
995cm- V A 



ping transport. The former corresponds to a delocalized 
picture, while the latter is to the localized one. Espe- 
cially, since the localization length in pentacene is es- 
timated to be about 10 molecules,»S this value could be 
directly utilized to estimate the decoherence time in pen- 
tacene. Within the present theory, we can also calculate 
the diffusion length within as the average distance be- 
tween initial and final site that before and after carrier 
hops. Obviously, when td is very large, our result should 
be the same with Troisi'sJ^ In Fig. |T1 it is found that, 
with td — 0.5 — Ips, the carrier could diffuse over about 
5 sites, that is, the spatial extent of the carrier is about 
10 times intermolecular distance, which is very close to 
the experimental prediction in pentacene 4 Hence, we can 
safely say that, Ips is the typical order of decoherence 
time in pentacene. Besides, with small td, the carrier is 
only localized within the next-nearest sites, such that, 
hopping transport should be the dominant for this case. 
Correspondingly, with temperature increasing, the large 
td curves decrease first and then increase slightly, while 
the small td curve does almost not change. This means, 
for the bandlike transport, the spatial extent of carri- 
ers should decrease with increasing temperature, while 
for the hopping transport, the carrier is mainly localized. 
The results for other parameters are similar, so we do not 
show here. 

The physical picture of these results is quite different 
from that in the traditional understanding of incoherent 
hopping, where the carrier residing in a trap would tun- 
nel to another one with similar energy. '^^ This mechanism 
works based upon the simple phonon structure and dilute 
impurities in inorganic semiconductors. In organic semi- 
conductors, however, due to the dynamic disorder each 
molecule might be treated as a trap and the phononic en- 
vironment is quite diverse. So the basic point for the dy- 
namics disorder model is, in our opinion, not the energy 
difference between molecules, but how long the carrier 
spreads. Fig. [S] shows a brief schematic for the present 
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FIG. 5: Schematic for the transition between hopping and 
bandlike transport. Under low temperature, the configura- 
tion of molecules is regular, while under high temperature, 
it is disordered. For small td, the carrier is mainly localized 
in one molecule, such that temperature acts as an assistance. 
For large td, the wavepacket of carrier is extended, and tem- 
perature produces disorders for the carrier's diffusion. 



theory. When td is small, the hopping is mainly between 
two neighbor molecules. In this case, due to the ther- 
mal motion of each molecule, the intermolecular distance 
eventually decreases when the temperature is sufficiently 
high, such that the hopping rate should be larger than 
that under low temperature. This is actually the temper- 
ature assisted hopping. While when td is large, that is, 
the carrier is much easier to spread out to several sites, 
the thermal motion of molecules behaves now as the dis- 
ordered obstacles. So in this case, temperature will act 
to produce more and more scattering to the carriers and 
then block the carrier's motion. 



D. Diffusion coefficient versus td 

In the present theory, the decoherence time td is the 
most important parameter. In Fig. [51 we show the de- 
pendence of diffusion coefficient on td- As expected, the 
diffusion coefficient increases following the increasing td, 
and when td > 2ps, diffusion coefficient tends to saturate 
to the result from Troisiii^ This is easy to understand, 
since when td is extremely small, it is quite hard for a 
carrier to hop out of the initial site. It is worth not- 
ing that, the estimated value of mobility in dynamics 
disorder model is slightly larger than that from exper- 
iment, say 3cm^s~^V^^ from theory and Icm-^s^^V^^ 
from experiment. The present theory based on the de- 
coherence might be the solution to this difference. Fur- 
thermore, the experiment on organic field effect transis- 
tor has shown that, when the drain voltage is large, the 
temperature coefficient is negative, and for the inverse 
case, it is positive.^ This could also be explained in the 
present theoretical framework. In case the drain volt- 
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age is large, the waiting time in each molecule should 
be small. Considering the ratio between and td that, 
the smaller the i^,, the larger the t^, the relationship of 
mobility with drain voltage could be easily understood. 

IV. CONCLUDING REMARKS 

Actually, the influence of phonons on charge transport 
in organic semiconductors has been studied extensively in 
the literature. However, the meaning of "incoherent hop- 
ping" in organic materials is still obscure in our opinion, 
since the widely used M-A formula is not such applica- 
ble for organic molecules. Meanwhile, as we show in this 
work, hopping and bandlike transport happens at differ- 
ent time scales. The first case refers to decoherence time 
of the same order with relaxation time of phonon so that 
phonon is able to transfer its energy to the carrier to as- 
sist it hopping, while the second case functions when td 



is large enough to ensure the time for scattering between 
carriers and phonons. In a real material, there are many 
ways to adapt these two conditions and make the type 
of transport transit from hopping to bandlike. Further- 
more, under low temperature, tunneling might become 
dominating, and carriers will search for a molecule out 
of the trap which has closest energy to tunnel to, i.e., 
Mott's variable range hopping mechanism. ■^^ This case 
has not been addressed in the present work. 

In the end, we would like to discuss more on the present 
theory. Firstly, the practical device parameters, such as 
electric held, have not been explicitly included, which 
could be easily considered within the present framework. 
For example, the electric held could be regarded as a 
phase factor in the hopping constant;^ but for the com- 
mon Ehrenfest method, this way loses its efficiency due 
to the Blocli oscillation. Especially, based on the study 
of electric field, we might discuss the validity of Einstein 
relation in these systems.'*'^ Secondly, it is straightfor- 
ward to replace the dynamic disorder model discussed 
in the present work, such as Holstein model, spin-boson 
model, etc. Especially, since the spin motion in organic 
materials is always realized to be coherent, it could be 
easily incorporated into the present theory, and the hot 
debated magnetic field effect could be studied^^ Thirdly, 
the theory could be applied to estimate the correspond- 
ing decoherence time in the specific material, since it is 
the only adjustable parameter. 
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